Cavity  Expansion  Analysis 
of  Brittle  Materials 


S.  Satapathy  and  S.  Bless 
Institute  for  Advanced  Technology 

The  University  of  Texas  at  Austin 


July  1995 


lAT.R  0075 


Approved  for  public  release;  distribution  unlimited. 


19960917  018 


Certification  of  Technical  Review 


Dr,  Thomas  Kiehne 


The  views,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the  author(s) 
and  should  not  be  construed  as  an  official  Department  of  the  Army  position,  policy,  or 
decision,  unless  so  designated  by  other  documentation. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  NO.  0704-0188 


Public  reporring  burden  tor  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  tor  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  anv  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services.  Directorate  for  Information  Operations  and  Reports,  1215  lefferson  Davis  Highway,  Suite  1204,  Arlington,  VA 
22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0180),  Washington,  DC  20503. 

1 .  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE 

3.  REPORT  TYPE  AND  DATES  COVERED 

July  1995 

Technical  Report 

4.  TITLE  AND  SUBTITLE 

Cavity  Expansion  Analysis  of  Brittle  Materials 

5.  FUNDING  NUMBERS 

Contract  #  DAAA21-93-C-0101 

6.  AUTHOR(S) 

S.  Satapathy  and  S.  Bless 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES} 

Institute  for  Advanced  Technology 

The  University  of  Texas  at  Austin 

4030-2  W.  Braker  Lane,  #200 

Austin,  TX  78759 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

lAT.R  0075 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  i 

U.S.  Army  Research  Laboratory 

ATTN:  AMSRL-WT-T 

Aberdeen  Proving  Ground,  MD  21005-5066 

10.  SPONSORING  /  MONITORING  AGENCY 

REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be  considered 
as  an  official  Department  of  the  Army  position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 

1  2a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

12b.  DISTRIBUTION  CODE 

A 

13.  ABSTRACT  (Maximum  200  words) 

In  this  report,  we  show  that  the  ‘target  resistance’  of  brittle  materials  can  be  calculated  accurately  using  spherical 
cavity  expansion  analysis  and  a  conventional  brittle  material  model.  The  stress  field  ahead  of  the  tip  of  the 
penetrator  is  assumed  to  have  spherical  symmetry.  The  brittle  material  is  modeled  as  an  elastic  material  which 
cracks  under  tension.  The  cracked  material  is  considered  to  be  pulverized  when  it  fails  in  compression,  which 
is  then  characterized  as  a  Mohr-Coulomb  material  with  pressure  dependent  shear  strength.  The  target  resistance 
value  found  from  this  analysis  compares  well  with  the  reported  experimental  values  for  AD995  alumina  (ALO3) 

14.  SUBIECT  TERMS 

target  resistance,  cavity  expansion  analysis,  brittle  materials,  ceramics, 
alumina. 

15.  NUMBER  OF  PAGES 

1  5 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT  . 

Unclassified 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

Unclassified 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

Unclassified 

20.  LIMITATION  OF  ABSTRACT 

UL 

NSN  7540-01-280-5500  Standard  Form  298  (Rev. 2-89) 


Prescribed  by  ANSI  Std.  239-18  298-102 


CAVITY  EXPANSION  ANALYSIS  OF  BRITTLE  MATERIALS 

by 

Sikhanda  Satapathy  and  Stephan  Bless 


1>0  Introduction 

Cavity  Expansion  Models  (CEM)  have  been  used  to  calculate  the  target  resistance  (Rt)  of 
materials  to  penetration.  Rt  is  considered  to  be  the  static  pressure  required  to  open  up  a 
cavity  in  the  target  [1,2],  This  approach  has  been  applied  successfully  for  metal  targets  first 
by  Bishop,  Hill,  and  Mott  [3]  and  more  recently  by  Forrestal,  et  al.  [4].  The  merit  of  this 
form  of  closed  form  analytical  solution  is  that  implications  of  material  models  can  be  directly 
determined  without  having  to  first  insert  them  into  wavecodes  and  numerically  simulate 
penetration  events.  This  kind  of  macroscopic  analysis  is  also  helpful  in  performing 
parametric  studies  which  can  identify  important  material  parameters  that  affect  penetration 
resistance. 

There  are  very  few  cavity  expansion  treatments  of  brittle  materials.  Forrestal,  et  al.  [5]  used 
an  elastic-cracked-plastic  material  model  to  calculate  Rt  for  AD85  alumina.  Wright,  et  al.  [6] 
used  an  elastic-plastic-hackle  material  model  to  calculate  Rt  for  polycarbonate.  Partom  used 
an  elastic-plastic  material  model  with  a  pressure  dependent  strength  in  the  plastic  region, 
which  saturates  at  high  pressure,  to  calculate  Rt  for  AD995  [7].  In  this  paper  we  demonstrate 
that  an  elastic-cracked-pul verized  material  model  yields  Rt  values  close  to  experimentally 
measured  values  for  AD995  alumina  [8].  We  do  not  assume  any  plastic  flow  in  the  target 
material.  It  is  found  that  the  target  resistance  depends  on  compressive  and  tensile  strengths 
and  the  pressure-shear  coefficient  of  the  pulverized  material. 


2.0  Response  Regions 

In  steady  state  penetration  (after  the  initial  shock  phase  is  over),  the  penetrator  erodes  and 
decelerates  due  to  the  stress  at  the  target-penetrator  interface.  The  stress  in  the  target  that 
resists  penetration  consists  of  two  parts:  1)  the  inertial  stress  required  to  accelerate  target 
material  out  of  the  path  of  the  penetrator;  and  2)  the  stress  required  to  deform  the  target 
around  the  cavity.  This  second  term  is  identified  with  the  target  resistance  Rt.  Partom  and 
Littlefield  [9]  pointed  out  that  there  is  also  a  dependence  of  Rt  on  velocity  in  materials  for 
which  the  strength  depends  on  pressure,  since  the  pressure  is  higher  in  high  velocity 
penetration.  That  coupling  is  ignored  in  this  treatment. 

The  stress  field  ahead  of  the  penetrator-target  interface  can  be  approximated  to  have 
spherical  symmetry.  We  model  the  alumina  ceramic  as  an  elastic  material  which  cracks 
under  tension  and  then  pulverizes  when  the  compressive  stress  exceeds  the  one  dimensional 
compressive  failure  strength.  The  implication  of  this  assumption  is  depicted  schematically  in 
Figure  1  which  shows  three  zones:  the  region  near  the  cavity  is  comminuted;  next  to  the 
comminuted  zone  there  is  a  “radially  cracked”  zone  (radial  cracks  are  formed  because  the 
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hoop  stress  exceeds  the  tensile  strength);  and  beyond  the  “radially  cracked”  zone  the  ater  i 
is  elastic.  This  phenomenology  is  consistent  with  Collombet  and  Tranchet’s  [10]  ob  vati 
of  explosive  cavities  in  alumina  (see  Figure  2);  they  find  “intergranular  cracks”  ^ar  t 
cavity,  next  to  which  “transgranular  cracks”  were  found.  The  “transgranular  cr'  ,s” 
radially  oriented  which  implies  that  for  spherical  symmetry,  the  hoop  stress  is  ze  In 
“intergranular  crack”  region,  the  cracks  are  randomly  oriented  and  are  inter  meets.  . 
Essentially,  the  material  in  this  region  is  fragmented  and  comminuted.  This  5  ^uence  of 
failed  regions  around  the  cavity  is  generally  consistent  with  measurements  aluminas 
penetrated  by  long  rods  [11].  Although  there  are  some  observations  th:  e  level  of 

communition  varies  in  the  pulverized  zone  [12],  we  model  the  entire  pulveri  aterial  as  a 

Mohr-Coulomb  material  with  a  pressure  dependent  shear  strength.  Th  .issure-shear 
coefficient  is  derived  from  the  uniaxial  stress-shear  coefficient  experimentally  obtained  by 
Clifton,  et  al.  [13]. 


Figure  1.  Response  regions  in  the  target 
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(a) 

Figure  2(a).  Isotropic  and  anisotropic  networks  of  cracks  [10].  2(b)  Transgranular 
cracks  in  the  soft  recovered  sample  [10]. 


View  A-A 


(b) 


3.0  Elastic  Region 

The  elastic  zone  extends  from  r  =  c  to  r  =  b  where  b  is  the  stress-free  outer  boundary  of  the 
target.  At  the  cracked-elastic  boundary  the  hoop  stress  equals  the  tensile  strength  of  the 
material.  We  neglect  inertial  terms,  since  Rt,  by  definition,  is  the  stress  required  to  open  a 
cavity  quasi-statically.  The  field  equations  can  be  written  as. 


Equilibrium  equation: 


dr  r 

du 


Strain-displacement  relation:  £r-~  ^^d  £  =  £  =  — 

dr  r 

Constitutive  relation  for  elastic  region:  e,  =  — [cT,.  -  2  vcr^] 

E 

WCT^- VO-,] 

Equations  (2),  (3),  and  (4)  can  be  combined  with  (1)  to  yield, 

d^u  2  du  2u 


A  z"*'— j - T  =  ® 

dr  r  dr  r 


(1) 

(2) 

(3) 

(4) 


(5) 


which  has  a  solution  of  the  form: 


u  =  c^r  +  -^ 
r 


(6) 
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To  evaluate  c,  and  we  need  two  boundary  conditions  which  are: 


Gr  =0 

^  r~b 


;  r 


7) 


and 


’e\r=c  - 


By  using  (2),  (3),  and  (4),  stresses  can  be  expressed  as: 


cr,  = 


'■  (1+v)(1-2v)L 
E 


.du  .  u 
(1-  V) —  +  2v- 
dr  r 


(^9  = 


du  u 

(1+ v)(1-2v)L'''^'^7J 

Thus,  using  (6)  to  (10),  the  constants  are  found  to  be: 

q‘^(i-2v) 


Cl  -  -y- 


,3  \ 


1  + 


2c' 


C2  =  — 


Gf(\-\-  v)b^ 
1,3  \ 


2E 


1  +  ■ 


2c' 


(8) 


(9) 

(10) 


(11) 


(12) 


Inserting  these  constants  in  (6),  (9),  and  (10)  yields  the  following  form  for  the  displacement 
and  stresses: 


u  = 


2E\ 


1  + 


2c' 


-  ,  a+v)b^ 
(l-2v)r  +  - - Y — 


(13) 


a.  = 


2c' 


„3 


(14) 


1  + 


1  + 


V 

2? 


2c' 


(15) 


Expressions  for  the  pressure  and  the  shear  stress  in  this  elastic  region  are  given  by: 
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(16) 


P  = 


1  +  — T 

2c^j 


<y. 


4 


1+ 


2c 


3 


y 


(17) 


4.0  Cracked  Region 

As  discussed  before,  the  material  will  crack  when  the  tensile  hoop  stress  exceeds  the  tensile 
strength.  Thus  the  elastic  region  will  be  bounded  inwardly  by  a  radially  cracked  region  (a  <  r 
<  c)  within  which  the  radial  cracks  are  assumed  to  be  spatially  distributed  in  various  sizes  in 
such  a  manner  that  the  hoop  stress  in  the  whole  region  is  zero.  The  inner  boundary  of  this 
region  is  defined  by  the  boundary  condition  that  the  material  is  pulverized.  Since  the  stress 
state  in  the  cracked  region  is  uniaxial,  pulverization  occurs  where  the  radial  stress  exceeds 
the  uniaxial  compressive  strength,  Y;  i.e.. 


(18) 

To  be  able  to  use  continuum  equations  in  the  cracked  region  (which  is  in  fact  discontinuous), 
we  resort  to  the  averaging  theorems  of  micro-mechanics  (see  Appendix  I)  which  let  us  apply 
the  continuum  equations  adopting  averaged  field  variables.  Hence  the  stress  and 
displacements  are  volume  averaged  quantities  allowing  for  local  variations.  Since  cr.  =0, 
the  equilibrium  equation  ( 1 )  reduces  to. 


dr 


which  can  be  integrated  with  (18)  to  yield. 


(19) 


(20) 


Also  since  the  radial  stress  has  to  be  continuous  at  the  elastic -cracked  boundary  r  =  c  then 
from  (14)  and  (20), 


1  + - T 

_ 2c^ 


3  A 


(21) 


To  evaluate  displacement  in  the  cracked  region,  Forrestal,  et  al.  [5]  have  used  the  equation  of 
state  with  bulk  modulus  of  the  intact  material.  But  in.  fact,  presence  of  radial  cracks  will 
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render  the  material  anisotropic.  Thus  we  can  not  use  the  isotropic  elastic  constants  of  the 
intact  material  to  characterize  the  displacement  of  this  anisotropic  cracked  region.  However, 
absence  of  transverse  cracks  ensures  no  degradation  of  elastic  moduli  in  the  radial  direction. 


'  Hence  using  (20)  with  the  stress-strain  relation  in  the  radial  direction,  (7,  =  Ee^ 
find. 


E  r 


A 


,22) 


The  constant  A  is  evaluated  by  using  continuity  of  displacement  at  r  =  c.  Thus  from  (13)  and 
(22),  we  obtain: 


u 


fj, 

\\-2v)c^  +{\+v)b^' 

2E 

2c') 

+ 


E  c  \r  ) 


The  pressure  and  shear  stress  in  this  region  can  be  written  as, 


P  = 


Ya^ 

2r^ 


(23) 


(24) 

(25) 


5.0  Comminuted  Region 

Subject  to  the  same  kind  of  assumptions  regarding  applying  continuum  equations  to  an 
inherently  discontinuous  region,  we  set  the  shear  stress  in  the  comminuted  region  (h  <  r  <  a) 
equal  to  the  “flow  stress,”  which  is  taken  as  pressure  dependent.  Thus, 


c  -cr^ 


•  =  m 


f  <7.  +  la  a 


(26) 


The  equilibrium  equation  (1)  reduces  to. 


^  +  2a-^  =  0 
dr  r 


where 


«:  = 


6m 


3  +  4m 


(27) 

(28) 


With  the  boundary  condition  (18),  (27)  can  be  solved  to  yield, 


6 


The  target  resistance,  Rt,  which  is  the  negative  of  the  radial  stress  at  the  cavity  boundarv 
r  =  h,  is  thus  found  as,  j  3^ 


(30) 


Finally,  the  pressure  and  the  shear  stress  in  the  comminuted  region  can  be  calculated  as, 


p  =  r 


r. 


Kr 


2a 


3-2a 


t  = 


(31) 

(32) 


6.0  Calculation  of  h/a 


From  (30),  we  note  that  to  solve  for  Rt,  we  need  to  evaluate  the  size  of  the  pulverized  zone 
relative  to  the  cavity  size.  Even  though  this  ratio  turns  out  to  be  a  constant,  independent  of 
geometry  for  an  infinite  target,  we  first  derive  its  value  for  a  finite  target.  Thus  far  we  have 
avoided  using  any  constitutive  relation  for  the  comminuted  region.  In  the  previous  section, 
we  integrated  the  equation  of  motion  with  the  assumption  that  the  shear  stress  varies  linearly 
with  pressure  in  this  region  (Mohr-Coulomb  criteria).  Even  if  we  assume  statistical 
homogeneity,  the  elastic  moduli  may  be  different  from  those  of  the  intact  material. 
Holmquist  [14]  notes  that  one  can  ignore  the  difference  between  the  elastic  moduli  of  intact 
and  pulverized  material.  Presence  of  bulking  due  to  comminution  complicates  the  problem 
further,  due  to  the  competition  between  dilatancy  and  pore  compaction.  Curran,  et  al.  [12] 
observed  that  dilatancy  occurs  only  after  the  large  compressive  stress  has  disappeared.  Thus 
if  we  neglect  both  dilatancy  and  compaction  due  to  pressure  (assuming  that  they  cancel  out  at 
least  to  a  first  order  approximation),  we  have  no  volume  change  in  the  comminuted  region. 
From  the  conservation  of  mass,  and  ignoring  higher  order  terms  in  u(a). 


i.f  ^  Y  _ 

3\aJ  a 


(33) 


Since  the  radial  displacement  has  to  be  continuous  at  r  =  a,  u(a)  is  given  by  (23).  Thus, 


l(h 

iU'  = 


2£’ 


3  \ 


1  + 


b 

2c^ 


(\-2v)c^  +  (\+v)b^ 


c^a 


+  -  1 


(34) 
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Equation  (34)  and  (21)  have  two  unknowns  c  and  a.  For  a  given  geometry  (h  and  b  kr 
these  two  equations  can  be  solved  for  “a”  and  the  cavity  expansion  pressure,  Rt,  t 
calculated. 


■^i), 

be 


7.0  InHiiite  Target 

For  an  infinite  target  (i.e.,  b  is  very  large),  the  RHS  of  (34)  is  a  function  of  only  c  vhich 
can  be  evaluated  from  (21).  Thus,  the  ratio  h/a  and  hence  Rt  is  determined  pureb  )m  the 
material  constants.  By  letting  b  -»  ©o,  (21),  (34)  and  (30)  reduce  to: 


c 

a 


\aj 


(1-v) 

V  T  V2 


3r, 


(l-v) 


^/2 


2a 

T 


(35) 


(36) 


(37) 


8.0  Material  Constants  for  Coors  AD995  Alumina 

We  calculate  the  elastic  constants  for  AD995  alumina  from  the  longitudinal  and  transverse 
wave  speeds  reported  by  Grady  [15]  as  E  =  373.14  GPa,  K  =  231.8  GPa.  The  quasi-static 
strength  parameters  as  published  by  Coors  are:  compressive  strength,  Y  =  2.62  GPa;  and 
tensile  strength  <t^  =  0.262  GPa.  Tracy  [16]  points  out  that  conventional  compression  tests 
tend  to  give  strength  measurements  that  are  misleadingly  low.  Using  a  dumbbell-shaped 
specimen  he  measured  a  compressive  strength  of  3.5  GPa  for  AD94  alumina.  As  a 
conservative  estimate,  compressive  strength  for  AD995  should  be  at  least  3.5  GPa.  Bar 
impact  tests  of  similar  alumina  gives  intermediate  results  [17,  18].  We  report  Rt  for 
compressive  strength  of  both  2.6  GPa  and  3.5  GPa.  Finally,  to  calculate  the  pressure-shear 
coefficient,  m,  we  make  use  of  Clifton’s  [13]  experimental  value  of  normal  stress  vs.  shear 
stress  coefficient  as  follows. 

Let  =  (38) 

solving  for  Gq,  we  obtain, 

Gg  =  G/l-2ju)  (39) 

and  hence  m  is  given  by: 
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(40) 


3/1 

{(J^  +  2Gq)/3  'i  —  AjX 

From  Clifton’s  [13]  measured  value  of  ii=  0.2,  we  obtain  m  =  0.273. 


9.0  Results 

We  plot  the  Rt  value  for  an  infinite  target  of  AD995  alumina  as  a  function  of  the  pressure 
shear  coefficient,  m,  in  Figure  3.  For  m  =  0.273  (from  Clifton’s  data),  we  find  that  Rt  =  7.7 
GPa  for  Y  =  2.62  GPa,  and  Rt  =  9.46  GPa  for  Y  =  3.5  GPa.  Subramanian  and  Bless  [8] 
conducted  penetration  experiments  on  AD995  targets  with  tungsten  penetrators  and  found 
that  Tate’s  [19]  solution  using  Rt=  8.5  GPa  passes  through  almost  all  the  error  bars  for  the 
experimental  values.  Thus,  our  closed  form  solution  has  resulted  in  Rt  values  that  tightly 
bracket  the  experimental  observation. 


Figure  3.  Variation  of  Rt  with  pressure-shear  coefGcient,  m. 


The  radial  stress  and  the  shear  stress  are  plotted  as  a  function  of  normalized  radius  (r/h)  in 
Figure  4.  It  can  be  seen  that  the  pulverized  zone  and  the  cracked  zone  extend  to  about  3.5 
and  9  times  the  cavity  size,  respectively.  The  stresses  fall  off  to  zero  asymptotically  beyond 
^^/h  >  15.  Thus  the  assumption  of  an  infinite  target  is  valid  if  the  target  to  penetrator  diameter 
ratio  is  at  least  15.  This  is  consistent  with  radial  confinement  effects  as  reported  by  Bless,  et 
al.  pO].  We  also  plot  in  Figure  5,  pressure  vs.  shear  stress  in  the  various  regions  which 
depict  the  path  taken  by  a  virgin  material  as  the  cavity  pressure  starts  to  build  up. 
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Figure  4.  Variation  of  stresses  with  normalized  radius,  r/h. 
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Figure  5.  Shear  stress  vs.  pressure  in  different  zones. 


10.0  Discussion 


This  relatively  simple  treatment  has  yielded  remarkably  good  correlation  with  measured  Rt 
values.  This  is  in  spite  of  several  weaknesses  in  the  approach,  namely:  a)  uncertainty  of  the 
measured  value  of  compressive  strength;  b)  accuracy  of  the  assumption  that  material  is 
uniformly  comminuted  in  the  pulverized  zone;  c)  uncertainty  of  measured  value  of  pressure- 
shear  coefficient  and  validity  of  using  this  value  (which  is  measured  in  alumina  powder)  for  a 
material  which  has  been  pulverized  from  AD995;  d)  inherent  assumption  of  CEM  that  Rt  is 
due  to  the  pressure  required  to  set  up  a  spherically  symmetric  stress  field  ahead  of  the 
penetrator;  e)  assumption  of  no  volume  change  in  the  pulverized  region;  and  f)  method  of 
calculating  Rt  from  experimental  data  by  arbitrarily  assuming  a  penetrator  strength  value  and 
assuming  the  Tate  formulation  [19]  correctly  accounts  for  inertial  stresses. 

We  feel  that  the  cavity  expansion  method  for  study  of  ceramic  penetration  has  considerable 
merit  and  should  be  pursued  further. 
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Appendix  I 


Frequently  we  have  to  deal  with  bodies  which  have  cracks,  voids,  inclusions,  etc.  All  real 
materials  have  this  kind  of  inhomogenity,  but  we  intend  to  neglect  these  local  variations  and 
apply  continuum  field  equations  using  remote  stress-strain  conditions.  Averaging  theorems 
let  us  use  the  continuum  equations  with  volume  averaged  quantities. 

Let  us  consider  a  representative  volume  element  (RVE)  where  external  traction  is 
prescribed  on  the  boundary  which  induces  the  local  stress  fields,  Oy.  Let  the  stress  field  as 
measured  on  the  boundary  be  CTy .  Now, 


j  Oydv  =  j  a■^5^^Jdv 

V 

=  J ^dv  ,  (since  in  absence  of  body  force,  cr,^  =  0) 

=  J  Gf^tij^Xjds,  (using  divergence  theorem) 
s 

=  \t^Xkds 

S 


=  J  aijUjX^ds 

s 


^GijV 

Hence, 


Similarly,  using  a  prescribed  displacement  field  on  the  boundary,  it  can  be  shown  that  the 
remote  strain  field,  ,  is  related  to  the  local  strain  field,  Sy,  as  follows: 

"  2y  I ^  "  F I 
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